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ABSTRACT 

We present the solution of the discrete super- Virasoro constraints to all orders 
of the genus expansion. Integrating over the fermionic variables we get a represen- 
tation of the partition function in terms of the one-matrix model. We also obtain 
the nonperturbative solution of the super- Virasoro constraints in the double scaling 
limit but do not find agreement between our flows and the known supersymmetric 
extensions of KdV. 
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1. Introduction. 

Matrix models provide us with powerful methods to evaluate the sum over ge- 
ometries in two dimensional quantum gravity. For c < 1 conformal field theories [1] 
coupled to gravity, n-point correlation functions of gravitationally dressed primary 
fields have been calculated within this formulation on world sheets with arbitrary 
topology, a result that is far from being obtained from the continuum Liouville 
approach [2] . It turns out that matrix models have a close connection to integrable 
systems [3], which is manifest through the appearence of the Korteweg-de Vries 
hierarchy in the double scaling limit [4] [5]. The recent discovery of Kontsevich ma- 
trix integrals has shown a natural connection to topological gravity [6] [7]. Despite 
of this important results, the generalization to the supersymmetric case is con- 
nected with difficulties. So for example, the role which the known supersymmetric 
extensions of KdV [8] play in this context remains obscure [9] [10] . For N = 1 
superconformal field theories coupled to two dimensional supergravity correlation 
functions for spherical topologies have been obtained in [11], using the continuum 
super-Liouville approach. The torus path integral and the spectrum of physical 
operators for c < 1 superconformal models coupled to 2D-supergravity have been 
considered in [9] [12]. However the more powerful counterpart in terms of super- 
matrix models has not yet been found [13]. To shed some light on this situation, 
the authors of [14] proposed an "eigenvalue model", which is described in terms 
of N even and iV-Grassmann odd variables. The basic idea for the construction 
of the model was the generalization of the Virasoro constraints known from the 
ordinary one matrix model [15] [16] to the N — 1 supersymmetric case. Gener- 
alizing Kazakov's loop equations [17] to N — 1 superloop equations , the model 
was evaluated for spherical topologies in [14], obtaining agreement with correlation 
functions [11] and critical exponents [18] [19] of (2,4m) superconformal field theo- 
ries coupled to world-sheet supergravity. Further analysis of the model carried out 
in [20] [21] showed, that the continuum super- Virasoro constraints are described in 
terms of a c = 1 theory, with a ^-twisted scalar and a Weyl Majorana fermion in 
the Ramond sector. The solution of the continuum superloop equations in the first 
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few orders of the genus expansion and a general heuristic argument showed [20], 
that there might be a close relation between the bosonic piece of the free energy 
and the free energy of the one matrix model. 

In this letter we analyze the solution of this model to all orders of the genus 
expansion. After a brief review about its construction (section two), we analyze 
the solution of the discrete super- Virasoro constraints in section three. Integrating 
over the fermionic variables, we get a representation of the partition function in 
terms of the one matrix model. In section four we analyze the solution of the 
super- Virasoro constraints in the double scaling limit and compare our flows with 
the known supersymmetric extensions of KdV [8]. 

2. Construction of the model. 

The bosonic one matrix model is defined as: 

ZB(9k,N)= [d N2 <S>exp[-^-trJ29k® k ] , A B = e"" (1) 

where $ is a Hermitian N x N matrix and \x is the bare cosmological constant. It 
satisfies the Virasoro constraints [15]: 

L n Z B = , n > -1 ; [L n , L m \ = (n - m)L n+m 
r A sf d " TTu 9 

n ~^h + h w 

which can be derived from (1) by implementing invariance under the shift $ — > 

$ + e $n+l. 

We obtain a relation to a free massless boson, if we introduce an infinite set of 
creation and annihilation operators related to the coupling constants 

N A B V2 d 
a- n = — j=ng n , n > ; a n = ^—t— , n > . (3) 

The Virasoro generators L n = \ ^ '■ a -k&k+n '■ are the modes of the energy 
momentum tensor T(p) = ^L n p~ n ~ 2 of the scalar field dX(p) = J2 a nP~ n ^ 1 = 
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dX + (p) + dX (p). From (1) and (3) we see that the potential of the model is 
related to the negative modes of the scalar field 



Jim \z B - x dX~Z B ~ V'(p) = V kg k p k -\ (4) 

N^oo iV z — ' 

fc>0 

As a generalization to the N — 1 supersymmetric case, we take a c"= 1 free mass- 
less superfield X(p, IT) = x(p) + 11^ (p) with energy momentum tensor T(p,U) = 
DXdX = ipdpX + II : (d p xd p x + dipip) : and the mode expansion dX{p) = 
Y^dnp^ 71 ^ 1 , ip(p) = ^2b r p~ r ~ 1 / 2 . The relation to the c = 1 super-Virasoro al- 
gebra: 



[L m , L n ] = (n - m)L n+m H <5 n+ m,o , (5) 



Tfl 

[L m , G r ] = ( — - r)G m +r , (6) 



{G r , G s } = 2L r+s + \(r 2 - \)S r+s ,o , (7) 
and the coupling constants of the model is introduced, if we define the modes as: 



A * 9 N n i 9 

«P = — ' a -p = -^—P9 P , P = 0,1,2,... 
N dg p A s 



In terms of the coupling constants the Virasoro generators then become: 
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00 9 00 n <9 A 2 n <9 2 

= V %fc^ + V( 7T + O^r^T + ^ V p 7j + 

^ ^fc+n ^ 2 <9£ n +r 2N 2 ^ dg k dg n ^ k 



^£<2 + '>5P^ n> -~ 1 - (10) 

In analogy to the one matrix model the potential is determined from the neg- 
ative modes of the superfield DV(p, II) ~ DX~: 

V(p,U) = J2(9kP k + m+i/2P k ) • (11) 

fc>0 

This suggests the introduction of N even-and N Grassmann-odd eigenvalues, 
obtaining for the partition function: 

N ( N 

Z s (g k ,t k+ i,N) = / l[dX i de i A(\,e)ex P J^Tt^i + t k+ ^i) 



i=l \ " k>0 i=l 



'12) 



The generalization of the usual Vandermonde determinant, A(A, 9), is deter- 
mined by implementing invariance of (12) under the super- Virasoro constraints (9) 
(10), which yields a differential equation for the measure: 



whose unique solution up to an irrelevant constant is A = Yli<j(^i~ — QiQj)- 
The model that we want to solve in the large N limit is thus: 



4 



z s = / n dx ^ - ^ - x£ s>* A * + • 

J i=l i<j \ S k>0 i=l J 

(14) 

In (14) N is taken to be even in order to get a Grassmann even partition 
function. 

3. Solution of the discrete model. 

We analyze the solution of the discrete superloop equations [14] [20] or equiv- 
alently the discrete super- Virasoro constraints through the representation (14) of 
the partition function. We first consider the case without fermionic couplings 
— 0. The partition function can be written as a function of the Pfaffian of 
the antisymmetric matrix My = A" 1 , where then A y - = Aj — Xj, i ^ j, and the 
usual Vandermonde A(A) = flio'C^i — -\?) [14] : 



2N „„ r 2N 



/IN 
Y[dX l A(X)Pf N (\r j 1 )exp(-—J2J29k^) • (15) 
i=i s k i=i 

Here the product of the Vandermonde and the Pfaffian is totally symmetric 
under exchange of two eigenvalues. We want to show, that the relation of the 
partition function (15) to the bosonic partition function of the one matrix model 
(1) is given by: 



Zs(9k,C k+ i = 0,27V) = ^^Z 2 B (9k,N) . (16) 



To prove (16) we show, that the following equivalence for the measures holds: 
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Pf N (^) - X >) = W ^ ^ J A 2 (/)A 2 ( J) , (17) 

i<3 I, J 

Here we have introduced the notation I = (A^, . . . , A^), i\ < . . . < in (same for 
J), and S ' means symmetric permutations between the elements of / and J. The 
equality (16) is a consequence of (17), since the partition function (15) with 2N 
eigenvalues will then factorize into ( 2 N ) products of two independent one matrix 
models with iV eigenvalues. 

Obviously (16) holds for the case of two eigenvalues, so that to illustrate the 
situation we start with the first nontrivial case of four eigenvalues, N = 2. We will 
consider the right and left hand side of (17) as polynomials Q(\), -P(A) respectively 
in A = A4: 



P(A) = (A-A 2 )(A-A3)Ai 2 Ai3-(A-A 1 )(A-A 3 )Ai 2 A 2 3+(A-Ai)(A-A 2 )Ai3A 2 3 (18) 

Q(X) = i[(A - A 3 ) 2 A? 2 + (A - A 2 ) 2 A? 3 + (A - Ai)^] . (19) 

From P(Xi) = Q(\i) for i = 1, 2, 3 it follows P(A) = Q(A), VA. 

We show that (17) is correct by induction over N. For iV > 2 we set A = A 2 at_|- 2 
so that -P(A), Q(X) become polynomials of degree 2N — 2 in A. It is enough to 
show that they coincide for A = X2N+1 since the symmetry under the interchange 
of Ai and Xj guarantees P(Xi) = Q(Xj) for i — 1, . . . , 2N and thus Q(X) = P(X) 
VN. If we consider both polynomials of (17) for N — > N + 1 and A = A 2 at + i we 
obtain the desidered result by applying the induction hypothesis. 



Q(A 2iV+ i) = E S J ' J A 2 0OAV) n ( A * - A 2iv+i) 2 (Ai - X 2N+1 ) 2 

I, J ii<...<ijv 
ji<...<jjv 
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2N i 

=pf 2N (K j i )u^- x ^ n (^-wi) 2 n (^-A 2 iv + i) 2 = p(A 2 iv + i) 

i<j h<...<i N ji<-<jN | 

(20) 

We conclude that the bosonic part of the supersymmetric (2iV x 2N) model is 
proportional to the square of the corresponding (N x N) bosonic one matrix model 
(16). With this result we can already find the form of the discrete string equation 
for arbitrary genus. 

Writing Z = e N Fs , equation (16) implies, that the fermionic independent 
piece of the supersymmetric free energy is related to the free energy of the 
one matrix model T& (up to an irrelevant additive constant): 



Ff = f • (21) 

Since the dependence of the free energy on go is trivial J|j = — 1/A, (21) implies 
a rescaling of the cosmological constant A$ = 2A#. This relation can also be 
obtained, if we evaluate the bosonic piece of the Lq constraint (10): 



A s = 2f;% fc 9 Ae (-A|^ . (22) 

1 

Up to a factor of two the cosmological constant satisfies the same string equa- 
tion as in the one matrix model. Including the fermionic couplings, the following 
equality for the second order in fermions holds: 



d 2 Z s (2N) 1 (2N\ f d 2 Z B (N) ^ dZ B {N)dZ B {N) \ 

di k+h di n+h -2^\n) {dg k+1 dg„ B[IV dg k+1 dg n ) ^ ~ U) " 

(23) 

The proof of this relation is straightforward, using the relation: 
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2N „ 2N 



a,(3=l J i=l 



I,J 



N 



(J)A 2 (J) £ (A£»A£ - At +1 A- ) - (fc - n) , (24) 



a,/3=l 



which can be shown using equation (17). Formula (23) is purely bosonic in the 
sense, that after derivation we set the fermionic couplings to zero. 

Equations (21) and (23) imply that the free energy until the second order in 
fermionic couplings is given by: 



Since the solution to the continuum super- Virasoro constraints will be analyzed 
in a moment, it will then become clear that (25) is the complete expansion of the 
free energy and no dependence on more than two fermionic couplings appears. 

Equation (25) gives a representation of the model (12) in terms of the one 
matrix model since the fermionic variables 9% have been integrated out. We think 
that further analysis of this interesting relation to the one matrix model for the 
discrete theory could lead to the geometrical interpretation of the model and the 
derivation of Feynman rules. We hope to report on this results elsewhere and turn 
now to the solution of the super- Virasoro constraints in the double scaling limit. 

4. Solution to the super- Virasoro constraints in the double scaling 
limit. 

The continuum super- Virasoro constraints [20] [21] are described by the super- 
energy momentum tensor of a c = 1 superconformal field theory 




(25) 
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T F = \d<p(z)1>(z) , T B = 1 - : dip{z)dip{z) : + 1 - : d^{z)^{z) : +^ , 

Z (27) 

where ip(z) is a bosonic scalar field with antiperiodic boundary conditions 
and ip(z) is a Weyl-Majorana Fermion in the Ramond sector. The modes of the 
fermionic part of the energy momentum tensor Tp = 5 X) G n+ iz~ n ~ 2 are described 
in terms of the coupling constants of the model for n > — 1: 



G 



faro, .Y^/i., 1 ^ 5 , r o 5 , 5 , 2f 52 



n+, 4k 2 1 Z_.v 2 y/c dr n+fc+1 2 9t„ ^ ^ fl+fc+1 Z^dt k dr n _ k \ 

(28) 

and the bosonic modes are obtained by the anticommutation relations (7). The 
constraints imposed on the partition function Zs(t k ,r k ) = exp(Fs) are: 



G n+ i_Z s (t k ,T k ) = , n>-\ . (29) 

The equations which describe the model (14) in the double scaling limit are 
determined by (28) (29): 



k 2 DF s = D- 1 u-2tDt , (30) 
t = -^2(2k + l)t k R k [u] , r = -J2nRk[u} , (31) 

k>l k>0 

where D = Jjj-, r = r + y = 0^ is the first fermionic scaling variable, t is 
the renormalized cosmological constant and u is the bosonic piece of the two point 
function of the puncture operator (ctqcto), satisfying the KdV-flow equation 
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du 



DR n+1 [u] 



(32) 



Rk[u] are the Gel'fand-Dikii polynomials defined through the recursion rela- 
tions: 



R = i, R x = u , R 2 = (3m 2 + k 2 D 2 u), . . . 

To show that (30) (31) is the solution to (29) it is enough to show that G_i 
and Gi are satisfied, since the algebra (5)-(7) guarantees for all the remaining 
constraints. That this is indeed the case, can be shown by simple calculations, 
where we have used well known relations of the one matrix model theory like (33). 
The equations (30) (31) coincide with the solution presented in [14] [20] for the 
first few orders of the genus expansion. The free energy (30) can be written as a 
function of the free energy of the one matrix model, which is very similar to the 
relation (26) of the discrete model and has an expansion at most bilinear in the 
fermionic couplings 



The form of the solution obtained (30) (31) implies, that the bosonic piece of 
the model is equivalent to an ordinary one matrix model, where the constraints 
act on the square root of the partition function L u \/Zq = 0, n > — 1 [16]. 
The fermionic dependence of the free energy determines, that correlation functions 
involving more than two fermionic operators vanish, since the fermionic scaling 
variable r is linear in fermionic couplings. These properties were known to hold 
for the first orders of the genus expansion [14] [20] and are in agreement with the 
previous analysis of the discrete model. 



DR n+1 [u] = (k 2 D 3 + 4uD + 2(Du))Rn[u] , 



(33) 




(34) 



k>i 
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The correlation functions in the double scaling limit which follow from (30) (31) 
as a function of the two point function of the puncture operator U = u — 2tD 2 t 
are given by: 



r)TJ <-> f)TJ f) 

^- = 2D(R k [u]Dr) , ^ = D R k+1 [u] - 2^-(tD 2 t) . (35) 

The scaling dimensions of the scaling operators (35) and the string suscepti- 
bility (31) are in agreement with the results of (2, 4m)-minimal superconformal 
models coupled to two dimensional supergravity [14] [18] [19]. 

From (35) we obtain for the first nontrivial bosonic flow, a set of two equations: 



|^ = D(3u 2 + k 2 D 2 u) , 1^- = QuDt + k 2 D 3 t , (36) 

where the first one is the ordinary KdV equation. Equations (36) are invariant 
under the global supersymmetric transformations 5u = eDr, 5r = eu, where e is a 
constant anticommuting parameter, as it is known to hold for the supersymmetric 
extensions of KdV of Manin-Radul or Mulase- Rabin [8]: 



= l^Wo + 3W d x W + \d 2 x W x W x , ^ = + \d x {W^W x ) 

(37) 

Here Wo is a bosonic variable, while W\ is fermionic. We can write the form of 
our generalized KdV-equation in terms of U, the perhaps more convenient variable: 



dU 

— = D(3U 2 + k 2 D 2 U - UDUtDt + 6k 2 DtD 3 t) , (38) 
at i 

8t 

= qUDt-12tDtD 2 t + k 2 D 3 t . (39) 
Equations (38) and (39) can independently be derived from the first orders of 
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the genus expansion, with the method proposed in [20]. From (36) or (38) (39) 
we have not been able to see any relation to the known supersymmetric extensions 
of the KdV-equation [8], so for example, the perhaps most natural identification 
Wq ~ U is not consistent with (37) nor the other extensions of KdV [8]. We 
may conclude that our correlation functions (35) describe a new supersymmetric 
extension of the KdV-hierarchy. 

The bosonic flows satisfy the recursion relations which can be obtained from 
the recursion relations of the i2fc[u]'s: 



du \ / K 2 D 2 + 2u + 2DuD- 1 n x /d 



u 



J±-J \ 2DTD- 1 + 2D- 1 Dt k 2 D 2 + 2u + 2D- 1 uD J 



(40) 



While for the fermionic flows we have simply: 



|^ = , 7 ^ = ^D 2 + 2u + 2D-'uD)^- (41) 
dr k dr k+1 dr k 



It would be important to study further properties of the differential equations 
which describe our correlation functions to see whether they share the standard 
properties of completely super-integrable systems. This will be left for work in the 
future. 
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